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Abstract 



c/5 ' In this work, in terms of suitable superpositions of equal-frequency Bessel beams, 

we develop a theoretical method to obtain nondiffractive beams in (weakly conduc- 
tive) absorbing media capable to resist the loss effects for long distances. 

q 
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>v 62.30.+d ; 43.60.+d ; 91.30.Fn ; 04.30.Nk ; 42.25.Bs ; 46.40.Cd ; 52.35.Lv . 
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^ : 1 Introduction 

o. . 

\Q ■ Over many years, the theory of nondiffracting beams and pulses has been developed, 

O _ 

generalized and experimentally verified in many fields, such as optics, microwaves and 

c*2 \ acoustics f 1-1 11. 

O 

When propagating in a non-absorbing medium, these waves maintain their spatial 



shape for long distances, i.e., they possess a large depth of field. 

However, the situation is not the same when dealing with absorbing media. In these 
cases, both the ordinary and the nondiffracting beams (and pulses) will suffer the same 
effect: an exponential attenuation along the propagation axis. 

Here, we are going to show that, through suitable superpositions of equal- frequency 
Bessel beams, it is possible to obtain nondiffracting beams in absorbing media capable to 
resist the loss effects, maintaining amplitude and spot size of their central core for long 
distances. 
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The method that will be developed here is a generalization of a previous one (also 
developed by us) O E3] , which was conceived for lossless media, where suitable super- 
positions of Bessel beams are used to construct stationary wave fields with arbitrary 
longitudinal intensity shape. Those new solutions were called Frozen Waves. 

Before continuing, some points are to be stressed: (a) the method developed here ad- 
dresses to poorly conductive media (a/eu « 1) in the optical frequency range; actually, 
we are considering materials where the penetration depth of the optical beams are typi- 
cally > l(r 4 m; (b) the energy absorption by the medium continues to occur normally, 
the difference being in that these new beams have an initial transverse field distribution, 
such to be able to reconstruct (even in the presence of absorption) their central cores for 
distances considerably longer than the penetration depths of ordinary (nondiffracting or 
diffracting) beams. 

In the next Section we shall extend the previous method fl2\ IT3] to absorbing me- 
dia. Section 3 is devoted to obtaining examples of these new beams in some interesting 
situations. 

2 The Mathematical Methodology 

The method that will be developed in this Section is based on Bessel beams superpo- 
sitions; so, it is appropriate first to recall some characteristics of such beams in absorbing 
media. 

In the same way as for lossless media, we construct a Bessel beam in the absorbing 
materials by superposing plane waves whose wave vectors lie on the surface of a cone 
with vertex angle 9. The refractive index of the medium can be written as n(u) = 
riR{uj) + ini(u), quantity being the real part of the complex refraction index and 
nj the imaginary one, responsible for the absorbing effects. With a plane wave, the 
penetration depth 5 for the frequency ujq is given by 5 = c/2ujoni. 

In this way, a zero-order Bessel beam in dissipative media can be written as 



2 



ip = J (k p p)exp(iPz)exp(—iujt) with f3 = n(u)uj cos 9/c = n R u cos9/c + injuj cos9/c = 
Pr + ifli] k p = n R u sin9/c + iniUJ sin#/c = k pR + ik pI , and so k 2 p = n 2 u 2 /c 2 — f3 2 . In 
this way ip = Jo((k pR + ik p i)p)exp(i(3 R z)exp(— iut)exp(— fliz), where (3r, k pR are the real 
parts of the longitudinal and transverse wave numbers, and k p j are the imaginary 
ones, while the absorption coefficient of a Bessel beam with an axicon angle 9 is given by 
ag = 2/3/ = 2njoj cos 6*/c, its penetration depth being 5 = 1/a = c/2un I cos9. 

Due to the fact that k p is complex, the amplitude of the Bessel function Jo(k p p) starts 
decreasing from p = till the transverse distance p = l/2k Pl , and afterwards it starts 
growing exponentially. This behavior is not physically acceptable, but one must remember 
that it occurs only because of the fact that an ideal Bessel beam needs an infinite aperture 
to be generated. However, in any real situation, when a Bessel beam is generated by finite 
apertures, that exponential growth in the transverse direction, starting after p = l/2k pn 
will not occur indefinitely, stopping at a given value of p. Actually, taking into account 
the size of the aperture generally used for generating optical Bessel beams, and the fact 
that we are considering media that are poor conductors, that exponential growth does 
not even happen, and the resulting Bessel beam only presents a decreasing intensity in 
the transverse direction. 

Let us now present our method. 

Consider an absorbing medium with the complex refraction index n{ui) = n R {uj) + 
ini(u)), and the following superposition of 2N + 1 Bessel beams with the same frequency 
u 

*(p,z,t) = £ A m J ((k pRm +zk pI Jp)e^^ z e- i ^ t e-^\ (1) 

m=-N 

where m are integer numbers, A m are constant (yet unknown) coefficients, /3 Rm and k pRm 
(Pi m and k p i m ) are the real parts (the imaginary parts) of the complex longitudinal and 
transverse wave numbers of the m-th Bessel beam in superposition (JTJ), the following 
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relations being satisfied 




= n 




(2) 



fa 



(3) 



01, 



Kir 



(4) 



where (3 m = fa m + ifa and k Pm = k pRm + ik pIm . 

Our goal is now to find out the values of the longitudinal wave numbers f3 m and the 
coefficients A m in order to reproduce approximately, inside the interval < z < L (on 
the axis p = 0), a freely chosen longitudinal intensity pattern that we call |-F(z)| 2 . 

The problem was already solved by us and the method developed for the particular 
case of lossless media |T21 E3], i.e., when rii = — > (3j m = 0. For those cases, it was 
shown that the choice f3 = Q + 2irm/L, with A m = J L F(z)exp(— i2irmz/L)/L dz can be 
used to provide approximately the desired longitudinal intensity pattern |F(2;)| 2 on the 
propagation axis, within the interval < z < L, and, at same time, to regulate the spot 
size of the resulting beam by means of the parameter Q, which can be also used to obtain 
large field depths (see details in Refs. fT^l H3]). 

However, when dealing with absorbing media, the procedure described in the last 
paragraph does not work, due to the presence of the functions exp(— (3i m z) in the super- 
position (^), since in this case that series does not became a Fourier series when p = 0. 

On attempting to overcome this limitation, let us write the real part of the longitudinal 
wave number, in superposition (JTJ), as 



fa 



Q + 



2irm 



(5) 



L 



with 



4 



< Q 



2irm lo 
—r- < n R — 
L c 



(6) 



where this inequality guarantees forward propagation only, with no evanescent waves. 
In this way the superposition ((!} can be written 

N 

9(p,z,t) = e-^'e"* £ A m J ((k pRm +ik pI Jp) e iS &' e^*-' , (7) 

m=—N 

where, by using Eqs.(JHJ), quantity j3i m is given by 

A„=^ = ^(Q + ^) (8) 



and k p m = k pR m + ik pi m is S iven b Y Eq.® . 

Now, let us examine the imaginary part of the longitudinal wave numbers. From 
Eq.(jHJ) the minimum and maximum values among the (3j m are (/3/) mm = (Q — 
2nN / L)ni /n R and (/3r) max = (Q + 27rN/L)nj/n R , the central one being given by 
/3 7 = (/3/) m=0 = Qnj/n R . With this in mind, let us evaluate the ratio 

Pi LQ 

Thus, when A << 1, there are no considerable differences among the various /3j m , 
since it holds /?j m ~ /3 7 for all m. And, in the same way, there are no considerable dif- 
ferences among the exponential attenuation factors, since exp(— (3i m z) ~ exp(— /3jz). So, 
when p = the series in the r.h.s. of Eq.0 can be approximately considered a truncated 
Fourier series multiplied by the function exp(— (3jz): and, therefore, superposition (JJJ) can 
be used to reproduce approximately the desired longitudinal intensity pattern |F(^)| 2 (on 
p = 0), within < z < L, when the coefficients A m are given by 
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(9) 



A m = - ( L F(z)e^ z e- i2 ^ z dz (10) 
L Jo 

the presence of the factor exp(/3jz) in the integrand being necessary to compensate for 
the factors exp(— /3i m z) in superposition 

Since we are adding together zero-order Bessel functions, we can expect a good field 
concentration around p = 0. 

In short, we have shown in this Section how one can get, in a weakly conductive ab- 
sorbing medium, a stationary wave-field with a good transverse concentration, and whose 
longitudinal intensity pattern (on p = 0) can approximately assume any desired shape 
|F(^)| 2 within the predetermined interval < z < L. The method is a generalization of a 
previous one O EH and consists in the superposition of Bessel beams in (JJJ), the real and 
imaginary parts of their longitudinal wave numbers being given by Eqs.(jS])and while 
their complex transverse wave numbers are given by Eq. ({%]). and, finally, the coefficients 
of the superposition are given by Eq.(fTUj). The method is justified, since 47iN/LQ « 1; 
happily enough, this condition is satisfied in a great number of situations. 

In the next Section we shall apply this method for obtaining new beams whose central 
core can resist the effects caused by diffraction and attenuation in absorbing media for 
long distances. 

3 Examples: Almost Undistorted Beams in Absorb- 
ing Media 

Let us consider a biological tissue as the absorbing medium. In the wavelength range 
280 — 300nm, proteins are the major absorbers, being the UV absorption similar for many 
types of tissues. 

A typical XeCl excimer laser (A = 308nm — > u = 6.09 x 10 15 Hz) has a penetration 
depth, in this case, of 5 cm; i.e., an absorption coefficient a = 20m _1 , and therefore 
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rij = 0.49 x 10~ 6 . The real part of the refraction index for this wavelength can be 
approximated by Ur = 1.5, and therefore n = + irij = 1.5 + i 0.49 x 10 -6 . 

A typical Bessel beam with ujq = 6.09 x 10 15 Hz and with an axicon angle 9 = 0.0141rad 
(so, with a transverse spot of radius 8.4 /im), when generated by an aperture, say, of radius 
R = 3.5 mm, can propagate in vacuum a distance (its field depth) equal to Z = Rj tan# = 
25 cm while resisting the diffraction effects. However, in the material medium considered 
here, the penetration depth of this Bessel beam would be only z p = 5 cm. 

Now, we can use the method of the previous Section to obtain, in the same medium 
and for the same wavelength, an almost undistorted beam capable of preserving its spot 
size and the intensity of its central core for a distance many times larger than the typical 
penetration depth of an ordinary beam (nondiffracting or not). 

With this purpose, let us suppose that, for this material medium, we want a beam 
(with u>o = 6.09 x 10 15 Hz) that maintains amplitude and spot size of its central core for 
a distance of 25 cm, i.e., a distance 5 times greater than the penetration depth of an 
ordinary beam with the same frequency. We can model this beam by choosing the desired 
longitudinal intensity pattern |F(z)| 2 (on p = 0), within < z < L, 



and by putting Z = 25 cm, with, for example, L = 33 cm. 

Now, the Bessel beam superposition can be used to reproduce approximately this 
intensity pattern, and to this purpose let us choose Q = 0.9999a;o/c for the /3^ m in Eq.©, 
and iV = 20 (which is allowed by inequality (JBJ))- 

Once we have chosen the values of Q, L and N, the values of the complex longitudinal 
and transverse Bessel beams wave numbers happen to be defined by relations (JEI) and 

Eventually, we can use Eq. fjlOJ) and find out the coefficients A m of the fundamental 
superposition 0, that defines the resulting stationary wave- field. 




elsewhere , 



1 for < z < Z 



(11) 
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Let us just note that the condition AnN/LQ << 1 is perfectly satisfied in this case. 

In Fig. 1(a) we compare the desired longitudinal intensity function |F(z)| 2 with the 
resulting stationary wave field, obtained from the Bessel beam superposition (JZj), and we 
can notice a good agreement between them. A better result could be reached by using 
a higher value of N, which in this specific example may assume, according to inequality 
©, a maximum value of 158. 

It is interesting to note that, at this distance (25 cm), an ordinary beam would have 
got its initial field-intensity attenuated 148 times. 




150 200 

z [mm] 



Figure 1: Comparison between the intensity of the desired longitudinal function F(z) 
(solid line) and that obtained from the Bessel beam superposition (JJJ) (dotted line). One 
can observe a good agreement between them. 



Figure 2(a) shows the 3D field-intensity of the resulting beam. One can see that the 
field possesses a good transverse localization (with a spot size smaller than 10/mi), it 
being capable of maintaining spot size and intensity of its central core till the desired 
distance. Figure 2(b) shows the same picture, but in an orthogonal projection. 
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Figure 2: (a) Three-dimensional field-intensity of the resulting beam, (b) The same 
figure in orthogonal projection. 

As we have said in the Introduction, the energy absorption by the medium contin- 
ues to occur normally; the difference is that these new beams have an initial transverse 
field distribution sophisticated enough to be able to reconstruct (even in the presence of 
absorption) their central cores till a certain distance. For a better visualization of this 
field- intensity distribution and of the energy flux, Fig. 3 (a) shows the resulting beam, in 
an orthogonal projection and in logarithmic scale. It is clear that the energy comes from 
the lateral regions, in order to reconstruct the central core of the beam. Figure 3(b) 
presents the initial field-intensity distribution on the aperture plane, at z — 0. We can 
note the complicated initial transverse field-intensity distribution needed for feeding the 
central core of the beam along the propagation axis. 

Obviously, the amount of energy necessary to construct these new beams is greater 
than that necessary to generate an ordinary beam in a non-absorbing medium. 
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Figure 3: (a) The resulting beam, in an orthogonal projection and in logaritmic scale, 
(b) Initial field-intensity distribution on the aperture plane, at z — 0. 

And it is also clear that there is a limitation on the depth of field of these new beams. 
For distances longer than 8 or 10 times the penetration depth of an ordinary beam, 
besides a great energy demand, we meet the fact that the field-intensity in the lateral 
regions would be even higher than that of the core, and the field would loose the usual 
characteristics of a beam (transverse field concentration). 

4 Conclusion 

In this paper we have shown that it is possible to construct almost-undistorted beams 
in absorbing media. They are got by suitable Bessel beams superpositions, and are capable 
to resist the loss effects, maintaining amplitude and spot size of their central core for long 
distances, when compared with the usual penetration depths of ordinary beams. 



10 



5 Acknowledgements 



The author is very grateful to Erasmo Recami, Hugo E. Hernandez-Figueroa and Claudio 
Conti for continuous discussions and kind collaboration. 

References 

[1] J. Durnin, J. J. Miceli and J. H. Eberly, "Diffraction-free beams," Physical Review 
Letters 58, 1499-1501 (1987). 

[2] I. M. Besieris, A. M. Shaarawi and R. W. Ziolkowski, "A bi-directional traveling 
plane wave representation of exact solutions of the scalar wave equation," J. Math. 
Phys. 30, 1254-1269 (1989). 

[3] J.-Y. Lu and J. F. Greenleaf , "Nondiffracting X-waves: Exact solutions to free-space 
scalar wave equation and their finite aperture realizations," IEEE Trans. Ultrason. 
Ferroelectr. Freq. Control 39, 19-31 (1992). 

[4] P. Saari and K. Reivelt, "Evidence of X-shaped propagation-invariant localized light 
waves," Phys. Rev. Lett. 79, 4135-4138 (1997). 

[5] E. Recami, "On localized X-shaped Superluminal solutions to Maxwell equations", 
Physica A 252, 586-610 (1998), and refs. therein. 

[6] D. Mugnai, A. Ranfagni and R. Ruggeri, "Observation of superluminal behaviors in 
wave propagation," Phys. Rev. Lett. 84, 4830-4833 (2000). This paper aroused some 
criticisms, to which the authors replied. 

[7] M. Zamboni-Rached, E. Recami and H. E. Hernandez-Figueroa, "New localized Su- 
perluminal solutions to the wave equations with finite total energies and arbitrary 
frequencies," European Physical Journal D 21, 217-228 (2002). 



11 



[8] For a review, see: E. Recami, M. Zamboni-Rached, K. Z. Nobrega, C. A. Dartora, and 
H. E. Hernandez-Figueroa, "On the localized superluminal solutions to the Maxwell 
equations," IEEE Journal of Selected Topics in Quantum Electronics 9, 59-73 (2003), 
and refs. therein. 

[9] S. Longhi, "Spatial-temporal Gauss-Laguerre waves in dispersive media", Physical 
Review E, 68, article no.066612 (2003). 

[10] M. A. Porras, S. Trillo, C. Conti, and P. Di Trapani, "Paraxial envelope X waves", 
Optics Letters, 28, 1090-1093 (2003). 

[11] C. Conti, and S. Trillo, "X waves generated at the second harmonic", Optics Letters, 
28, 1251-1253 (2003). 

[12] M. Zamboni-Rached, "Stationary optical wave fields with arbitrary longitudinal 
shape by superposing equal frequency Bessel beams: Frozen Waves," Optics Express 
12, 4001-4006 (2004). 

[13] M. Zamboni-Rached, E. Recami and H. E. Hernandez-Figueroa "Theory of Frozen 
Waves: Modelling the Shape of Stationary Wave Fields," J. Opt. Soc. Am., A, (to 
be published). 



12 



